Out-of-plane impurities induced the deviation from the monotonia d-wave 
superconducting gap in cuprate superconductors 
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The electronic structure of cuprate superconductors is studied within the kinetic energy driven 
superconducting mechanism in the presence of out-of-plane impurities. With increasing impurity 
concentration, although both superconducting coherence peaks around the nodal and antinodal re- 
gions are suppressed, the position of the leading-edge mid-point of the electron spectrum around the 
nodal region remains at the same position, whereas around the antinodal region it is shifted towards 
higher binding energies, this leads to a strong deviation from the monotonic d-wave superconducting 
gap in the out-of-plane impurity-controlled cuprate superconductors. 
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I. INTRODUCTION 



The superconducting gap is a fundamental property of 
superconductors^, and the nature of its anisotropy has 
played a crucial role in the testing of the microscopic 
theory of superconductivity in cuprate superconductors^. 
Experimentally, by virtue of systematic measurements^, 
particularly using the angle-resolved photoemission spec- 
troscopy (ARPES) technique^, the d-wave nature of the 
superconducting gap has been well established by now. 
In particular, this d-wave superconducting symmetry re- 
mains one of the cornerstones of our understanding of 
the physics in cuprate superconductors^iii^i^ii. The early 
ARPES measurements on the cuprate superconductor 
Bi2Sr2CaCu208+*^ showed that in the real space the gap 
function and the pairing force have a range of one lattice 
spacing, and then the superconducting gap function is of 
the monotonic d-wave form Ak = A[cosfc2: — cosky]/2. 
Later, the ARPES measurements on the cuprate super- 
conductor Bi2Sr2CaCu208-i-,^ indicated that the super- 
conducting gap significantly deviates from this mono- 
tonic d-wave form. Furthermore, it was argued that 
this deviation should be attributed to an increase of 
the electron correlation, which may increase the inten- 
sity of the higher order of the harmonic component in 
the monotonic d-wave gap function^. However, recent 
ARPES measurementsiSiii on the cuprate superconduc- 
tors (Bi,Pb)2(Sr,La)2Cu06+5 and Bi2Sri.6in.o.4Cu06+5 
(in-Bi2201) with Ln=La, Nd, and Gd showed that a 
much stronger deviation from the monotonic d-wave su- 
perconducting gap form is unlikely to be due to the strong 
correlation effect. 

The cuprate superconductors have a layered structure 
consisting of the two-dimensional Cu02 layers separated 
by insulating layers"^^'^^. The single common feature is 
the presence of the Cu02 plane^iH, and it seems evi- 
dent that the unusual behaviors of cuprate superconduc- 
tors are dominated by this Cu02 plane^. It has been 
well established that the Cu^"*" ions exhibit an antifer- 
romagnetic long-range order in the parent compounds 
of cuprate superconductors, and superconductivity oc- 
curs when the antiferromagnetic long-range order state 
is suppressed by doped charge carriers^^. Since these 
doped charge carriers are induced by the replacement 
of ions by those with different valences or the addition 



of excess oxygens in the block layer, therefore in prin- 
ciple, all cuprate superconductors have naturally impu- 
rities (or disorder). However, for the cuprate super- 
conductors (Bi,Pb)2(Sr,La)2Cu06+5 and Lrt-Bi2201, the 
mismatch in the ionic radius between Bi and Pb or Sr 
and Ln causes the out-of-plane impurities^, where the 
concentration of the out-of-plane impurities is controlled 
by varying the radius of the Pb or Ln ions, and then 
the superconducting transition temperature Tc is found 
to be decreasing with increasing impurity concentration. 
These cuprate superconductors (Bi,Pb)2(Sr,La)2Cu06+5 
and Lrt-Bi2201 are often referred to as the out-of-plane 
impurity-controlled cuprate superconductors. Recently, 
the electronic structure of the out-of-plane impurity- 
controlled cuprate superconductors and the related su- 
perconducting gap function have been investigated ex- 
perimentally by using ARPES ^°d^ . It was shown that 
although the effect of the out-of-plane impurity scatter- 
ing around the antinodal region is much stronger than 
that around the nodal region, both superconducting co- 
herence peaks around the nodal and antinodal regions are 
suppressed. Furthermore, the magnitude of the deviation 
from the monotonic d-wave superconducting gap form in- 
creases with increasing impurity concentrationjiSiii. The 
appearance of the strong deviation from the monotonic 
d-wave superconducting gap form observed recently in 
the out-of-plane impurity-controlled cuprate supercon- 
ductors (Bi,Pb)2(Sr,La)2Cu06+5 and Ln-Bi2201 is the 
most remarkable effec t^'^d^ , however, its full understand- 
ing is still a challenging issue. To the best of our knowl- 
edge, this strong deviation from the monotonic d-wave 
superconducting gap form in the out-of-plane impurity- 
controlled cuprate superconductors has not been treated 
starting from a microscopic superconducting theory yet. 

In the absence of out-of-plane impurity scattering, 
the electronic structure of cuprate superconductors 
in the superconducting state has been discusse d^^d^ 
within the framework of the kinetic energy driven 
superconductivity^^, where the superconducting gap 
function has a monotonic d-wave form, and the main fea- 
tures of the ARPES experiments^ on cuprate supercon- 
ductors have been reproduced. In this paper, we study 
the electronic structure of the out-of-plane impurity- 
controlled cuprate superconductors in the superconduct- 
ing state and the related superconducting gap function 
along with this line. We employ the t-J model by con- 
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sidering the out-of-plane impurity scattering, and then 
show exphcitly that the strong deviation from the mono- 
tonic d-wave superconducting gap form occurs due to 
the presence of the impurity scattering. Although both 
sharp superconducting coherence peaks around the nodal 
and the antinodal regions are suppressed, the effect of 
the impurity scattering is stronger in the antinodal re- 
gion than that in the nodal region. Our results also show 
that the electronic structure of the out-of-plane impurity- 
controlled cuprate superconductors in the superconduct- 
ing state can be understood within the framework of the 
kinetic energy driven superconducting mechanism with 
the out-of-plane impurity scattering taken into account. 

This paper is organized as follows. In Sec. II we 
present the basic formalism of the electronic structure 
calculation in the presence of the out-of-plane impurities. 
Within this theoretical framework, we discuss the elec- 
tronic structure of the out-of-plane impurity-controlled 
cuprate superconductors in the superconducting state 
and the related superconducting gap function in Sec. 
Ill, where we show that the well pronounced deviation 
from the monotonic d-wave superconducting gap form is 
mainly caused by the out-of-plane impurity scattering. 
Finally, we give a summary in Sec. IV. 



II. FORMALISM 

It has been shown that the essential physics of 
cuprate superconductors is properly accounted by the 
two-dimensional t-J model on a square lattice^, 

H ^ ^ C/g. Ci+f,a + ^ C}^ Ci+fa + M ^ C}^ C'iff 

ifjcr if a ia 

+ jX]S,-S,+^, (1) 

ifj 

acting on the Hilbert subspace with no doubly occupied 
site, i.e., J^a ^L^i'y ^ 1' where f) = ±x, ±y, t = ±x±y, 
Cj^ (Cicr) is the creation (annihilation) operator of an 
electron with spin a, Si — {Sf , Sf , S^) are spin opera- 
tors, and /i is the chemical potential. To deal with the 
constraint of no double occupancy in analytical calcula- 
tions, the charge-spin separation fermion-spin theory— 
has been developed, where the constrained electron op- 
erators Cii and Cii are decoupled as Ci| = hl^S^ and 

Cii — hl^S^ , respectively, here the spinful fcrmion op- 
erator hi^ = e~***''/ii describes the charge degree of 
freedom together with some effects of spin configura- 
tion rearrangements due to the presence of the doped 
charge carrier itself, while the spin operator Si describes 
the spin degree of freedom, then the electron on-site lo- 
cal constraint for the single occupancy, J^a^ia^ic = 
S+h^^hl^S- + SrL,^hlS+ = hMiStS- + S~S+) = 

1 — hlhi < 1, is satisfied in analytical calculations. In 
particular, it has been shown that under the decoupling 
scheme, this charge-spin separation fermion-spin repre- 
sentation is a natural representation of the constrained 
electron defined in the Hilbert subspace without double 
electron occupancy^^. Furthermore, these charge carrier 



and spin are gauge invariant-^, and in this sense they are 
real and can be interpreted as physical excitations^. In 
this charge-spin separation fermion-spin representation, 
the t-J model (1) can be expressed as, 



H 



if] 
ir 

^J■ E h]a^i<^ + Jeff E • Si+fj, (2) 



with Jeff = (1 - S)'^J, and S = {h\„hi„) = be- 
ing the charge carrier doping concentration. This Jcs 
is similar to that obtained in Gutzwiller approach^. As 
an important consequence, the kinetic energy term in 
the t-J model has been transferred as the interaction be- 
tween charge carriers and spins, which reflects that even 
the kinetic energy term in the t-J Hamiltonian has a 
strong Coulombic contribution due to the restriction of 
no double occupancy of a given site. This interaction 
from the kinetic energy term is quite strong, and it has 
been showr^ii in terms of the Eliashberg's strong coupling 
theory^ that in the case without an antiferromagnetic 
long-range order, this interaction can induce a charge car- 
rier pairing state (then the electron Cooper pairing state) 
with d-wave symmetry by exchanging spin excitations in 
the higher power of the charge carrier doping concentra- 
tion 6. In this case, the electron Cooper pairs originat- 
ing from the charge carrier pairing state are due to the 
charge-spin recombination, and their condensation re- 
veals the d-wave superconducting ground-state. Further- 
more, this d-wave superconducting state is controlled by 
both the superconducting gap function and the quasipar- 
ticle coherence, which leads to the fact that the maximal 
superconducting transition temperature occurs around 
the optimal doping, and then decreases in both under- 
doped and overdoped regimea^i. Moreover, it has been 
shown^'"^'^^ that this superconducting state is the conven- 
tional Bardeen-Cooper-Schrieffer (BCS) likeii^ with the 
d-wave symmetry, so that the basic BCS formalism with 
the d-wave superconducting gap function is still valid in 
quantitatively reproducing all main low energy features 
of the ARPES experimental measurements on cuprate su- 
perconductors, although the pairing mechanism is driven 
by the kinetic energy by exchanging spin excitations, and 
other exotic magnetic scatterings, is beyond the BCS for- 
malism. Following previous discussionsi^ii^iii, the full 
charge carrier Green's function in the superconducting 
state with a monotonic d-wave gap function can be ob- 
tained in the Nambu representation as2^. 



5(k,w) 



^hF- 



i / tj + ^k A,,z(k) 

UJTq + A/tz(k)Tl -f ^kT3 



(3) 



where tq is the unit matrix, ti and ra are the Pauli matri- 
ces, the renormalized charge carrier excitation spectrum 
= ZhF£,k, with the mean- field charge carrier excitation 
spectrum = Ztxij-k — Zt'x2l\s_ — f-, the spin correla- 
tion functions xi = {S:t S^+fi) ^-nd X2 = {StSr^~), 7k = 
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(1/^) T,f, e'""^, 7k = (1/^) Ef e'*"'^, ^ is the number of 
the nearest neighbor or next nearest neighbor sites, the 
renormahzed charge carrier monotonic d-wave pair gap 
function AhzO^) = ZhF^hO^), where the effective charge 
carrier monotonic d-wave pair gap function A/i(k) = 

^hj\^^ with 7^'*' = (cosfcj: — cosfcj,)/2, and the charge 



carrier quasiparticle spectrum Ehu = YCk+ I ^/iz(k) p. 
The charge carrier quasiparticle coherent weight Z^f and 
effective charge carrier gap parameter Ah are determined 
by the foUowing two equationai^ii^iii, 



F«(k,p,p') (k,P,p') 



(2), 



k,p,p' 



Ehk ujpUjp, \ (iVpi -ajp^ - El^^ [uop, + ujpY ~ El^. 



(4a) 



pWp/ I (Wp - UJp. 



+ 



p,p 
i^f (P,P') 



^hp-p'+ko)^ ('^p - t^p' + £^hp-p'+ko)^ 



(p,p') 



+ 



{ujp + UJp' — Ehp-p'+ko)^ {cOp + uip' + Ehp-p'+ko 



(4b) 



respectively, where ko = [7r,0], Ak = Ztji^ — Zt'j[^, 
Bp = 2Ai(Ai7p - A2) - A2(2xi7p ~ X2), Ai = 2ZJ,s, 
A2 = 4Zcl)^t', A, = exl + Xi/2, A2 = xl + eXi/2, 
e = 1 + 2t0i/Jeff, the charge carrier's particle-hole 
parameters 4>i = {hl^hi+f,^) and 02 = {h\^hi+fa), 
the spin correlation functions xf = {SfSf_^f-^) and 

= (^f^f+f), Fi'''(k,p,p') = (^p, - c.p)[nB(c.p) - 
nB{u)p>)][l - 2nFiEhk)] + Ehk[nB{uJp')nB{-u;p) + 

nB(wp)nB(-cjp/)], Ff^(k,p,p') = (wp- + 
Wp)[n_B(-cjp/) - n_B(a;p)][l - 2n_F(£;,,k)] + 



Ehk[nB{i^p')nB{t^p)+nB{-UJp')nB{-UJp)], (p,p') = 
nF{Ehp^p,+k^)[nBiuJp') - UBiuJp)] - ns(cjp)n_B(-cjp'), 

-Fi^^(P,P') = nFiEhp-p>+i,o)[nB{oJp) 

"-b('^p')] - ns(t^p')^iB(-t^p), F2^^^(p,p') = 
nF{Ehp^p'+k^)[nB{uJp') - nB{-u}p)] + nB{uJp)nBiuJp'), 

-Fi'^^(P,P') = nFiEhp-p'+]^a)[nB{~ujp') - ns(wp)] + 
riB(— tjp)nB(— cjp/), ns(wp) and npiEhia) are the boson 
and fermion distribution functions, respectively, and the 
mean-field spin excitation spectrum, 



/ ,2 - \2 



A4 - aexi7p - j (1 - e7p) + Ms - ^"Xf - aXi7p ) (e " 7p) 



1 



77 Ms - 7raX2 



2Z 
+ A1A2 



X27p - ^X2 j 7p 



1 



axf(l ~ e7p)7p + 2"(Xi7p - C'3)(e - 7p) + a7p(C3 - 6X2 7p) - -^aeiCs - X27p) 



, (5) 



where A3 = aCi + {l-a)/{2Z), A4 = aCl + [l-a)/{AZ), 
A^ = aC2 + (1 — a)/{2Z), and the spin corre- 
lation functions C, = (1/^') E^,,' (^ir)*?- ^,), - 

= (l/^)Ef(^if,^r+f): and C| = 
(V^)Ef('5f+^^f+f)- In order to satisfy the sum 
rule of the correlation function (S^S^) — 1/2 in the 
case without the antiferromagnetic long-range order, an 
important decoupling parameter a has been introduced 
in the above calculatio n^^i^^''^^ , which can be regarded 
as the vertex correction. These two equations (4a) 
and (4b) must be solved simultaneously with other 
self-consistent equations, then all order parameters, the 
decoupling parameter a, and the chemical potential fi 
are determined by the self-consistent calculatio n^^i^^i^^ . 
In this sense, the calculations in this kinetic energy 
driven superconductivity scheme are controllable with- 



out using any adjustable parameters. We emphasize 
that the Green's function (3) is obtained within the 
kinetic energy driven superconducting mechanism, 
although the similar phenomenological expression has 
been used to discuss the impurity effect in cuprate 
superconductors2ii2^. 

With the charge carrier BCS formalism (3) under the 
kinetic energy driven superconducting mechanism, we 
can now introduce the effect of impurity scatterers into 
the electronic structure. In the presence of impurities, 
the unperturbed charge carrier Green's function in Eq. 
(3) is dressed by impurity scattering^^, 

g/(k,w)-i = g(k,a;)-i - Z,;^S(k,^), (6) 

with the self-energy function E(k, w) = I]ct(k, ti;)rQ,. 
In this case, the charge carrier Green's function in Eq. 
(6) can be explicitly rewritten as. 
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gi{k,uj) = 2^5/Q(k,w)Ta = ZhF 



[UJ - So(k,Cj)]To + [A,,^(k) + I]l(k,C^)]Tl + [a + S3(k,0^)]T3 

[c^ - So(k,c^)]2 - [Ck + S3(k,c^)]2 - [Ahz{k) + Si(k,c^)]2 



(7) 



Based on this Green's function (7), have discussed 
the effect of the extended impurity scatterers on the 
quasiparticle transport of cuprate superconductors in the 
superconducting state within the nodal approximation 
of the quasiparticle excitations and scattering processes, 
where the main effect on the quasiparticle transport 
comes from the extended impurity forward (or diagonal) 
scatterers, and therefore the component of the self-energy 
function Si(k, w) has been neglected, while the compo- 
nents of I]o(k, cj) and I]3(k, a;) have been treated within 
the framework of the T-matrix approximation. However, 
it has been demonstrated that the superconducting tran- 
sition temperature is considerably affected by the out- 
of-plane impurity scattering in spite of a relatively weak 
increase of the residual resistivity^. This reflects the 
fact that the superconducting pairing is very sensitive 
to the out-of-plane impurity scattering, and then the ef- 
fect of the out-of-plane impurity scattering is always ac- 
companied by breaking of the superconducting pairs. In 
this case, the out-of-plane impurities can be described 
as the elastic off- diagonal scatterers or pairing impurity 
scatterers. In particular, the modulation of the out-of- 
plane impurity scattering potential observed in scanning 
tunneling microscopy experiments^ has a characteristic 
wavelength of a few lattice spacings, this may arise be- 
cause the impurities give rise to an atomic-scale modula- 
tion of the charge carrier pairing potential which causes 
larger, coherence length size fluctuations in the out-of- 
plane impurity scattering potential^. Furthermore, the 
crude effect of the order parameter modulations on the 
quasiparticle scattering by allowing the order parameter 
to be modulated on the four bonds around the impurity 
has been estimated^- by adding the off-diagonal scatter- 
ing potential. 



obtaine d^^'^^ as V"(q) — piV{q)p{q), where 



V = 



Y,[Vik) + V{k')]n 

k,k' 

}-Vo y^}{coskx - cosky) + (cosfc^ - cosfc^)]ri, (8) 



k,k' 



to the phenomenological d-wave BCS Hamiltonian, then 
it was shown that the scattering rate is largest at the 
antinode. 

The exact form of the out-of-plane impurity scattering 
potential is very important for a better understanding 
of the electronic structure of the out-of-plane impurity- 
controlled cuprate superconductors. In the following 
discussions, we determine the form of the out-of-plane 
impurity scattering potential in terms of the calcula- 
tion of Dyson's equation. The potential which scatters 
the electron is taken as summation of impurity poten- 
tials V — J2i ~ where the summation is over 
all impurity sites /, and then its Fourier transform is 



p(q) 
Pi(q) 



(9) 
(10) 



are the charge carrier density in the Nambu representa- 
tion and the impurity density, respectively. In the cal- 
culation of the self-energy function induced by the im- 
purity scattering, usually it is assumed that the impuri- 
ties are randomly located and that there is no correla- 
tion between their positions^*'^^. In this case, the self- 
energy function can be obtained as S(k, uj) — &^'> (k, lu) + 
I]^2)(k, w) within the Born approximation, with the cor- 
responding first-order and second-order self-energy func- 
tions are evaluated a o^^'^^ , 

E(i)(k,c^) = p.^(5k'=onk') =P^^(0), (11a) 

k' 

f:(2)(k,c^) = p, 5k'+k"=o"^"(k')g/(k + k',c^)^(k") 

k',k" 

= p,^V^(k')57(k + k',c.)F(-k'), (lib) 



where pi is the impurity concentration. As we have 
mentioned above, the out-of-plane impurities are the off- 
diagonal scatterers. Although their scattering has a very 
weak effect on the residual resistivity for cuprate super- 
conductors, a heavy effect on the d-wave SC state is ob- 
served experimentally^*. With these considerations, we 
introduce the following out-of-plane impurity scattering 
potential, 

V = Y^ V{k')Ti = T^o ^[cosfc^ - cosfc;]Ti. (12) 



In this case, V^(0) ^ Fo[cos(0) - cos(O)] = (then 
Si(k,w) = 0), and S(k,w) = f;(2)(k,tj). This form of 
the out-of-plane impurity scattering potential in Eq. (12) 
is very similar to that in Eq. (8) used in Ref. |25| . and 
the scattering rate is also largest at the antinode. This 
is indeed confirmed by the quantitative characteristics 
presented in the following section. With the help of the 
impurity scattering potential in Eq. (12), the compo- 
nents of the charge carrier self-energy function E(k, w) 
are obtained explicitly as. 



I]o(k,c.) = p,l^|l^(k')p5/o(k' + k,u;) 

k' 

= A4El^(k'-k)p5,o(k',a;), (13a 
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E3(k,c.) = -p,l^|y(k')|V3(k' + k,c.) 

k' 

= -Pr^J2\^(^' -^^\^9i3ik',u;), (13b) 

k'+k 

I]i(k,c.) = p,l^|^(k')|V(k' + k,c.) 

k' 
k' 

In the charge-spin separation fermion-spin theorjsi^, 
the electron diagonal and ofF-diagonal Green's functions 



are the convolutions of the spin Green's functio n^^'^^i^' 
D~^{p,uj) = (w^ — LUp)/Bp and the charge car- 
rier diagonal and off-diagonal Green's functions in 
Eq. (7), respectively. These convolutions reflect the 
charge-spin recombination'^'^. Following the previous 
discussionsi^ii^iii, we can obtain the electron diagonal 
and off-diagonal Green's functions in the present case. 
Then the electron spectral function from the electron di- 
agonal Green's function is found explicitly as, 



MK^) = -T;^i^{[nBi^^p) + nFiujp-uj)]Ah{p-k,ujp-uj) 



N ^ 2wp 
p ^ 



[nB(-Wp) -I- n_F(-Wp - u;)]A/i(p - k, -cjp - w)}, (14) 
I 



where Ah is the charge carrier spectral function, which tained from Eq. (7) as, 
can be expressed as Ah = — 2Imgj*"(k, w), with 5*" ob- 



dia,. N 7 ^-So(k,^)+$k + S3(k,^) 

^ ' ^ '''^[L.-Eo(k,a;)]2-Kk + E3(k,c.)]2-[A;,^(k) + Si(k,a;)]2- ^ ' 

I 



III. ELECTRONIC STRUCTURE OF THE 
OUT-OF-PLANE IMPURITY-CONTROLLED 
CUPRATE SUPERCONDUCTORS 



Experimentally, it has been shown that the average of 
the next-nearest neighbor hopping t' is not appreciably 
affected by the out-of-plane impurities^J^ . In this case, 
the commonly used parameters in this paper are cho- 
sen as tj J = 2.5 and t' jt = 0.3. We are now ready 
to discuss the electronic structure of the out-of-plane 
impurity-controlled cuprate superconductors and the re- 
lated superconducting gap. In cuprate superconductors, 
the information revealed by ARPES experiments'' has 
shown that around the nodal [7r/2,7r/2] and antinodal 
[tt, 0] points of the Brillouin zone contain the essentials of 
the whole low energy quasiparticle excitation spectrum. 
In this case, we have performed a calculation for the elec- 
tron spectral function yl(k, w) in Eq. (14) at both nodal 
and antinodal points. The results at (a) the nodal point 
and (b) the antinodal point with the impurity concen- 
tration Pi — 0.001 (solid line), pi = 0.002 (dashed line), 
and Pi = 0.003 (dotted line) under the impurity scatter- 
ing potential with Vb = 50 J for the charge carrier doping 
concentration (5 = 0.15 are plotted in Fig. 1. For compar- 
ison, the corresponding ARPES experimental results^^ 
for the out-of-plane impurity-controlled cuprate super- 
conductors L7i-Bi2201 in the superconducting state are 



also presented in Fig. 1 (inset). Our results show that 
the quasiparticle peak is strongly dependent on the im- 
purity concentration, and the peaks at both nodal and 
antinodal points are suppressed due to the presence of 
impurity scattering. At the nodal point, there is a sharp 
superconducting quasiparticle peak near the Fermi en- 
ergy, however, although the peak at the high impurity 
concentration is dramatically reduced compared to that 
at the low impurity concentration, the position of the 
leading-edge mid-point of the electron spectral function 
remains almost unchanged. In particular, the position of 
the leading-edge mid-point of the electron spectral func- 
tion reaches the Fermi level, indicating that there is no 
superconducting gap. On the other hand, the spectral in- 
tensity from the Fermi energy down to ~ —1.1 J decreases 
as the impurity concentration increases at the antinodal 
point, this is the same case as that at the nodal point. 
However, the position of the leading-edge mid-point of 
the electron spectral function is shifted towards higher 
binding energies with increasing impurity concentration, 
this is in contrast with the behavior observed at the nodal 
point, and indicates the presence of the superconducting 
gap. The present results also show that the effect of the 
out-of-plane impurity scattering is stronger at the antin- 
odal point than at the nodal one, in qualitative agreement 
with the experimental resultsi^iii. 

The behavior of the electron spectrum in Fig. 1 in- 
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FIG. 1: The electron spectral function at (a) the nodal point 
and (b) the antinodal point with pi = 0.001 (solid line), pi — 
0.002 (dashed line), and pi = 0.003 (dotted line) for Vo = 50 J 
at 5 = 0.15. Inset: the corresponding experimental results 
taken from Ref. 
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(solid line) for Vo = 50J at 5 — 0.15. Inset: the corresponding 
experimental results taken from Ref. [lol . 



dicates an enhancement of the superconducting gap in 
the antinodal region by the impurity scattering. To show 
this point clearly, we have calculated the electron spectral 
function A(k, lj) along the direction [7r,0] [7r/2,7r/2], 
and then employed the shift of the leading-edge mid-point 
as a measure of the magnitude of the superconducting 
gap at each momentum just as it has been done in the 
experimentsiSiiii. The results for the extracted supercon- 
ducting gap as a function of the Fermi surface angle 0, de- 
fined in the inset, with the impurity concentration pi = 
(dashed line) and pi = 0.001 (solid line) under the im- 
purity scattering potential with Vq = 50 J for the charge 
carrier doping concentration (5 = 0.15 are plotted in Fig. 
2 in comparison with the corresponding ARPES exper- 
imental results for the out-of-plane impurity-controlled 
cuprate superconductor (Bi,Pb)2(Sr,La)2Cu06+5 in the 



"1 — I — I — r 
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■■■ Nd-Bi2201 
■ La-Bi2201 




0.2 0.4 0.6 0.8 
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FIG. 3: The superconducting gap as a function of [cosfc^, — 
cosfcj,]/2 with pi — 0001 (solid line), p,; = 0.002 (dashed line), 
and pi = 0.003 (dotted line) for Vo = 50J at 5 = 0.15. Inset: 
the corresponding experimental results taken from Ref. [Til . 



superconducting stateiS (inset). It is clearly shown that 
the superconducting gap A increases with the Fermi sur- 
face angle decreasing from 45° (node) to 0° (antinode). 
Although the superconducting gap in the presence of the 
impurity scattering is basically consistent with the d- 
wave symmetry, it is obvious that there is a strong devi- 
ation from the monotonic d-wave form around the antin- 
odal region. In particular, this strong deviation is mainly 
caused by a remarkable enhancement of the supercon- 
ducting gap value around the antinodal region, in qual- 
itative agreement with the experimental resultsiSiii. In 
other words, the superconducting gap around the antin- 
odal region is strongly enhanced by the impurity scatter- 
ing, whereas around the nodal region its value remains 
the same. As a consequence, the well pronounced devi- 
ation from the monotonic d-wave superconducting gap 
form in the out-of-plane impurity-controlled cuprate su- 
perconductors is mainly caused by the effect of the out- 
of-plane impurity scattering. This is also the reason why 
the superconducting gap function for very high quality 
samples of the cuprate superconductor Lai_a;Sra:Cu04 
has a monotonic d-wave form^. 

For a better understanding of the impurity concen- 
tration dependence of the deviation from the monotonic 
d-wave superconducting gap function, we have made a 
series of calculations for the superconducting gap at dif- 
ferent impurity concentration levels, and the results of 
the superconducting gap as a function of the monotonic 
d-wave function [cosfca; — cosfcy]/2 with the impurity con- 
centration Pi = 0001 (solid line), pi — 0.002 (dashed 
line), and pi = 0.003 (dotted line) under the impurity 
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scattering potential with Vq = 50 J for the charge carrier 
doping concentration S — 0.15 are plotted in Fig. 3 in 
comparison with the corresponding ARPES experimental 
results for the out-of-plane impurity-controlled cuprate 
superconductors Ln-Bi2201^^ (inset). Obviously, our re- 
sults show that the magnitude of the deviation from the 
monotonic d-wave superconducting gap form around the 
antinodal region increases with increasing impurity con- 
centration, in qualitative agreement with the experimen- 
tal resultsiSiii. This strong out-of-plane impurity effect 
in the antinodal region is also consistent with scanning 
tunneling spectroscopy results^, where the average of 
the superconducting gap size, which corresponds to the 
antinodal superconducting gap in the ARPES spectra, 
increases with increasing impurity concentration. 

Within the framework of the kinetic energy driven su- 
perconducting mechanisniii in the presence of the out-of- 
plane impurities, our present results show that the out- 
of-plane impurity scattering potential (12) in which the 
impurities modulate the pair interaction locally give qual- 
itative agreement with respect to the main features ob- 
served in the ARPES measurements on the out-of-plane 
impurity-controlled cuprate superconductors in the su- 
perconducting state. Although this out-of-plane impu- 
rity effect in cuprate superconductors can also be dis- 
cussed starting directly from a phenomenological d-wave 
BCS formalism^i^, in this paper we are primarily in- 
terested in exploring the general notion of the effects of 
the out-of-plane impurity scatterers in the kinetic energy 
driven cuprate superconductors in the superconducting 
state. The qualitative agreement between the present 
theoretical results and ARPES experimental data also 
indicates that the presence of the out-of-plane impurities 
has a crucial impact on the electronic structure of cuprate 
superconductors. On the other hand, we emphasize that 
the quasiparticle scattering rate in the antinodal region 
is strongly increased by the impurity scattering potential 
(12), while the nodal quasiparticles are very weakly scat- 
tered by the impurity scattering potential (12), this is 
why the superconducting transition temperature is con- 
siderably affected by the out-of-plane impurity scatter- 
ing in spite of a relatively weak increase of the residual 
resistivitjii^, since the transport properties are mainly 
governed by the quasiparticles in the nodal region. 



IV. SUMMARY 

In conclusion, we have shown very clearly in this paper 
that if the out-of-plane impurity scattering is taken into 
account within the framework of the kinetic energy driven 
d-wave superconductivityil, the quasiparticle spectrum 
of the t-J model calculated based on the off-diagonal im- 
purity scattering potential (12) per se can correctly repro- 
duce some main features found in the ARPES measure- 
ments on the out-of-plane impurity-controlled cuprate 
superconductors in the superconducting state^iii. In the 
presence of the out-of-plane impurities, although both 
sharp superconducting coherence peaks around the nodal 
and antinodal regions are suppressed, the effect of the 
impurity scattering is stronger in the antinodal region 
than that in the nodal region, this leads to a strong de- 
viation from the monotonic d-wave superconducting gap 
form in the out-of-plane impurity-controlled cuprate su- 
perconductors. 

Finally, we have noted that within a phenomenological 
BCS approach, the electron spectral properties of the un- 
derdoped cuprates as resulting from a momentum depen- 
dent pseudogap in the normal state have been discussed^, 
where a normal state pseudogap function deviating from 
the monotonic d-wave pseudogap form has been used to 
fit the ARPES experimental data in the normal state. It 
has been show n^^i'^^ that there are some subtle differences 
for different families of underdoped cuprates in the nor- 
mal state, and therefore it is possible that the pseudogap 
in the normal state is effected by the impurity scattering 
as well. 
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